Crimes_2001_to_present <- read_csv("Crimes_-_2001_to_present.csv", 
     col_types = cols(Latitude = col_character(), 
        Location = col_character(), Longitude = col_character(), 
        `X Coordinate` = col_character(), 
        `Y Coordinate` = col_character()))
Crimes_2001_to_present
Crimes <- Crimes_2001_to_present
Crimes08 <- Crimes %>% filter(Year=="2008")
Crimes08
Crimes08Type <- Crimes08 %>% group_by(`Primary Type`, `Location Description`) %>% summarise(count = n()) %>% rename(PriType=`Primary Type`, LocDes = `Location Description`)
Crimes08Type
Graph08 <- ggplot(data=Crimes08Type,aes(x=reorder(PriType,count),y=count ,fill=LocDes))+geom_bar(stat='identity',position='stack', width=.9)+theme(axis.text.x = element_text(angle = 45, hjust = 1))+scale_color_grey()
ggplotly(Graph08, height=1000, width=3000)
EachYear <- Crimes %>% filter(Year!=2020) %>% group_by(Year,`Primary Type`) %>% summarise(count=n()) %>% arrange(`Primary Type`) %>% rename(PriType=`Primary Type`)
EachYear
EachYearGraph <- ggplot(data=EachYear,aes(x=Year,y=count,fill=PriType))+geom_bar(stat='identity',position='stack', width=.9)+theme(axis.text.x = element_text(angle = 45, hjust = 1))+scale_color_grey()
ggplotly(EachYearGraph, height=1000, width=1500)
options(noaakey = "dZHCRTVsEdwFQEzgovApkrhRWtwHWjjJ")
ncdc(datasetid = 'GHCNDMS', stationid = 'GHCND:USW00014819', startdate = '2010-01-01',
   enddate = '2019-12-31')
$meta
$meta$totalCount
[1] 1200

$meta$pageCount
[1] 25

$meta$offset
[1] 1


$data

attr(,"class")
[1] "ncdc_data"
dat <- ghcnd(stationid = "USW00014819")
using cached file: /Users/kairui/Library/Caches/R/noaa_ghcnd/USW00014819.dly
date created (size, mb): 2020-04-24 21:35:18 (1.589)
dattemp <- dat %>%
  filter(element == "TMAX", year >= 2000)%>% select(VALUE1,year,month)%>%group_by(VALUE1=VALUE1/10)%>%group_by(AvgTemp=VALUE1)%>%select(year, month, AvgTemp)
dattemp
Crimes08 <- Crimes08 %>% mutate(Date = lubridate:: dmy_hms(Date))%>%arrange(Date)
 257218 failed to parse.
Crimes08
Crimes08 <- Crimes08 %>% group_by(`Primary Type`,month=floor_date(Date, "month")) %>% summarise(count = n()) %>% rename(PriType=`Primary Type`)
Crimes08
Crimes08Theft <- Crimes08 %>% filter(PriType=="THEFT")
Crimes08Theft
ggplot( data = dattemp, aes(x = AvgTemp)) + geom_density(adjust = 0.4) + labs(title = "")

LS0tCnRpdGxlOiAiRmlsZSIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKYGBge3IgaW5jbHVkZT1GQUxTRX0KbGlicmFyeShEYXRhQ29tcHV0aW5nKQpsaWJyYXJ5KHRpZHl2ZXJzZSkKbGlicmFyeShyZWFkcikKbGlicmFyeShtb3NhaWMpCmxpYnJhcnkobHVicmlkYXRlKQpsaWJyYXJ5KHBsb3RseSkKYGBgCgoKYGBge3J9CkNyaW1lc18yMDAxX3RvX3ByZXNlbnQgPC0gcmVhZF9jc3YoIkNyaW1lc18tXzIwMDFfdG9fcHJlc2VudC5jc3YiLCAKICAgICBjb2xfdHlwZXMgPSBjb2xzKExhdGl0dWRlID0gY29sX2NoYXJhY3RlcigpLCAKICAgICAgICBMb2NhdGlvbiA9IGNvbF9jaGFyYWN0ZXIoKSwgTG9uZ2l0dWRlID0gY29sX2NoYXJhY3RlcigpLCAKICAgICAgICBgWCBDb29yZGluYXRlYCA9IGNvbF9jaGFyYWN0ZXIoKSwgCiAgICAgICAgYFkgQ29vcmRpbmF0ZWAgPSBjb2xfY2hhcmFjdGVyKCkpKQpgYGAKCgpgYGB7cn0KQ3JpbWVzXzIwMDFfdG9fcHJlc2VudApgYGAKCmBgYHtyfQpDcmltZXMgPC0gQ3JpbWVzXzIwMDFfdG9fcHJlc2VudApDcmltZXMwOCA8LSBDcmltZXMgJT4lIGZpbHRlcihZZWFyPT0iMjAwOCIpCkNyaW1lczA4CmBgYAoKYGBge3J9CkNyaW1lczA4VHlwZSA8LSBDcmltZXMwOCAlPiUgZ3JvdXBfYnkoYFByaW1hcnkgVHlwZWAsIGBMb2NhdGlvbiBEZXNjcmlwdGlvbmApICU+JSBzdW1tYXJpc2UoY291bnQgPSBuKCkpICU+JSByZW5hbWUoUHJpVHlwZT1gUHJpbWFyeSBUeXBlYCwgTG9jRGVzID0gYExvY2F0aW9uIERlc2NyaXB0aW9uYCkKQ3JpbWVzMDhUeXBlCmBgYAoKYGBge3J9CkdyYXBoMDggPC0gZ2dwbG90KGRhdGE9Q3JpbWVzMDhUeXBlLGFlcyh4PXJlb3JkZXIoUHJpVHlwZSxjb3VudCkseT1jb3VudCAsZmlsbD1Mb2NEZXMpKStnZW9tX2JhcihzdGF0PSdpZGVudGl0eScscG9zaXRpb249J3N0YWNrJywgd2lkdGg9LjkpK3RoZW1lKGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KGFuZ2xlID0gNDUsIGhqdXN0ID0gMSkpK3NjYWxlX2NvbG9yX2dyZXkoKQpnZ3Bsb3RseShHcmFwaDA4LCBoZWlnaHQ9MTAwMCwgd2lkdGg9MzAwMCkKYGBgCgpgYGB7cn0KRWFjaFllYXIgPC0gQ3JpbWVzICU+JSBmaWx0ZXIoWWVhciE9MjAyMCkgJT4lIGdyb3VwX2J5KFllYXIsYFByaW1hcnkgVHlwZWApICU+JSBzdW1tYXJpc2UoY291bnQ9bigpKSAlPiUgYXJyYW5nZShgUHJpbWFyeSBUeXBlYCkgJT4lIHJlbmFtZShQcmlUeXBlPWBQcmltYXJ5IFR5cGVgKQpFYWNoWWVhcgpgYGAKCmBgYHtyfQpFYWNoWWVhckdyYXBoIDwtIGdncGxvdChkYXRhPUVhY2hZZWFyLGFlcyh4PVllYXIseT1jb3VudCxmaWxsPVByaVR5cGUpKStnZW9tX2JhcihzdGF0PSdpZGVudGl0eScscG9zaXRpb249J3N0YWNrJywgd2lkdGg9LjkpK3RoZW1lKGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KGFuZ2xlID0gNDUsIGhqdXN0ID0gMSkpK3NjYWxlX2NvbG9yX2dyZXkoKQpnZ3Bsb3RseShFYWNoWWVhckdyYXBoLCBoZWlnaHQ9MTAwMCwgd2lkdGg9MTUwMCkKYGBgCgpgYGB7ciBpbmNsdWRlPUZBTFNFfQpkZXZ0b29sczo6aW5zdGFsbF9naXRodWIoInJvcGVuc2NpL3Jub2FhIikKYGBgCgpgYGB7ciBpbmNsdWRlPUZBTFNFfQpsaWJyYXJ5KCJybm9hYSIpCmBgYApgYGB7cn0Kb3B0aW9ucyhub2Fha2V5ID0gImRaSENSVFZzRWR3RlFFemdvdkFwa3JoUld0d0hXampKIikKbmNkYyhkYXRhc2V0aWQgPSAnR0hDTkRNUycsIHN0YXRpb25pZCA9ICdHSENORDpVU1cwMDAxNDgxOScsIHN0YXJ0ZGF0ZSA9ICcyMDEwLTAxLTAxJywKICAgZW5kZGF0ZSA9ICcyMDE5LTEyLTMxJykKYGBgCgpgYGB7cn0KZGF0IDwtIGdoY25kKHN0YXRpb25pZCA9ICJVU1cwMDAxNDgxOSIpCmRhdHRlbXAgPC0gZGF0ICU+JQogIGZpbHRlcihlbGVtZW50ID09ICJUTUFYIiwgeWVhciA+PSAyMDAwKSU+JSBzZWxlY3QoVkFMVUUxLHllYXIsbW9udGgpJT4lZ3JvdXBfYnkoVkFMVUUxPVZBTFVFMS8xMCklPiVncm91cF9ieShBdmdUZW1wPVZBTFVFMSklPiVzZWxlY3QoeWVhciwgbW9udGgsIEF2Z1RlbXApCmRhdHRlbXAKYGBgCmBgYHtyfQpDcmltZXMwOCA8LSBDcmltZXMwOCAlPiUgbXV0YXRlKERhdGUgPSBsdWJyaWRhdGU6OiBkbXlfaG1zKERhdGUpKSU+JWFycmFuZ2UoRGF0ZSkKQ3JpbWVzMDgKYGBgCgpgYGB7cn0KQ3JpbWVzMDggPC0gQ3JpbWVzMDggJT4lIGdyb3VwX2J5KGBQcmltYXJ5IFR5cGVgLG1vbnRoPWZsb29yX2RhdGUoRGF0ZSwgIm1vbnRoIikpICU+JSBzdW1tYXJpc2UoY291bnQgPSBuKCkpICU+JSByZW5hbWUoUHJpVHlwZT1gUHJpbWFyeSBUeXBlYCkKQ3JpbWVzMDgKYGBgCmBgYHtyfQpDcmltZXMwOFRoZWZ0IDwtIENyaW1lczA4ICU+JSBmaWx0ZXIoUHJpVHlwZT09IlRIRUZUIikKQ3JpbWVzMDhUaGVmdApgYGAKCgpgYGB7cn0KZ2dwbG90KCBkYXRhID0gZGF0dGVtcCwgYWVzKHggPSBBdmdUZW1wKSkgKyBnZW9tX2RlbnNpdHkoYWRqdXN0ID0gMC40KSArIGxhYnModGl0bGUgPSAiIikKYGBgCgoKCg==